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Spati a 1 -frequency characteristics (SFC) and the scattering function 
of a layer which is uniform on the horizontal with absolutely black bottom, 
or an isolated layer, are the basis for a mathematical model which 
describes the horizontal heterogeneities in a light field with regard for 
radiation polarization in a three-dimensional planar atmosphere, delimited 
by a heterogeneous surface with diffuse reflection. The perturbation 
method was used to obtain vector transfer equations which correspond to 
the linear and nonlinear systems of polarization radiation transfer. The 
SFC of the nonlinear system satisfy the parametric set of one-dimensional 
boundary value tasks for the vector transfer equation, and are expressed 
through the SFC of linear approximation. As a consequence of the devel- 
oped theory, formulas were obtained for analytical calculation of albedo 
in solving the task of dissemination of polarization radiation in the 
planetary atmosphere with uniform Lambert bottom. 

Key words: vector transfer equation, polarization radiation, three- 

dimensional planar layer, heterogeneous albedo, spatial -frequency charac- 
teristics, scattering functions, fundamental solution, linear and nonlinear 
system. 
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POLARIZATION RADIATION IN THE PLANETARY ATMOSPHERE DELIMITED 
BY A HETEROGENEOUS DIFFUSELY REFLECTING SURFACE 

By S. A. Strelkov and T. A. Sushkevich 


Introduction /3* 

Study of the polarization characteristics of planetary atmospheric 
radiation has been underway already for over 100 years, if we even begin 
to count from the founding work [1], whose author for the first time 
proposed the four parameters which bear his name to describe the partially 
polarized Stokes 1 i ght-vector . Calculation of polarization radiation in 
repeatedly scattering heterogeneous, but even isotropic media with non- 
Rayleigh law of radiation interaction with substance, until now has pre- 
sented great computational difficulties, although a lot of attention of 
many specialists of different countries using the latest achievements in 
the field of modern numerical and analytical methods, and also the most 
powerful computers has been focused on these problems in the last 30 years. 
At the same time, a trend has been observed towards expanding the sphere 
of application of the polarization radiation in studies by specialists 
working in different fields, including astrophysics, atmospheric and space 
optics, hydrooptics, radio physics, ecology (environmental pollution), 
remote sensing of the earth's surface, ocean, seas, reservoirs, etc., in 
relation to studying natural resources. Naturally it is necessary to 
develop methods which allow calculation of the effect of heterogeneous 
underlying surfaces on the formation of the polarized radiation field of 
the medium. Since solution to the transfer equation with regard for 
polarization in itself is a cumbersome task, it is necessary to develop a M 
method for calculating the diffuse reflection from the heterogeneous 
underlying surface that does not result in additional complications and 
considerable increase in the volume of the computations. 

This publication suggests a mathematical model which corresponds to 
these requirements. It is based on one-dimensional vector kinetic equa- 
tions for spatial-frequency characteristics (SFC) of a system of polarized 
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radiation transfer , i nvari ant in relation to the constant component and 
horizontal variations in the underlying surface albedo and the conditions 
of layer illumination. The method of integrated Fourier transform was 
used to construct the fundamental solutions, or the functions of the 
influence [2] of a three-dimensional vector transfer equation with per- 
turbed boundary condition. In the theory of vision in turbid media, and 
in the tasks of Fourier optics, the functions of influence are usually 
called scattering functions [3, 4], The Stokes vector is presented in 
the form of an integrated packet of fundamental solutions or SFC with 
horizontal variations in the albedo of the underlying surface by the 
method of summing the perturbations. These results are the development 
and generalization for the vector task of the results presented in publica- 
tion [5, 6]. 

The formula for calculating the albedo of a uniform Lambert surface, 
depolarizing the radiation in the atmosphere follows from a constructed 
solution as a particular case. This correlation permits a significant 
reduction in the volume of computations, since calculation of the Stokes 
vector is required to solve the task of calculating the Stokes vector in 
the atmosphere-reflecting underlying surface system only for an isolated 
heterogeneous atmosphere. This formula was derived for the first time and 
is a generalization for the vector case of the formula for a scalar task 
of the transfer theory previously obtained by V. V. Sobolev [7, 8], Van 
de Holst [12], Z. G. Yanovitskiy [9], I. N. Minin [10], V. V. Ivanov [11] 
and the authors of publication [13]. 

In order to study the link between polarization radiation reflected 
by the planet, the characteristics of diffuse reflection of its surface, 
an expression was found for direct optical transfer operator of the atmos- 
phere, generalization of the vector task for the result presented in pub- 
1 ication [13]. 

1. Mathematical Statement of the Task 

t 

0 ' . 

An external parallel stream of radiation it S^t falls on the upper 

boundary of a planar layer unlimited in horizontal coordinates and finite 
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on the vertical (0<z<H, H-geometric thickness of the layer) at angle 
mO/K&iju* t0 the perpendicular to the layer with azimuth , £-■( 

The vector t describes the condition of polarization of the incident 
light streams; £ = { 1,0, 0,0} corresponds to nonpolarized, natural light. 
tiS x is the power of the stream, \ is the index of the wavelength of mono- 
chromatic radiation. Direction of dissemination of radiation t is 
described by spherical coordinates: [o,iJ K 

zenith angle which is counted from the direction of the inner perpendic- 
ular to the upper layer boundary z = 0. It coincides with the axis z, 
and --azimuth counted from the positive direction of the x- 

axis. We adopt as /*0, the azimuth of the external monodirectional 
stream , i/.* 0 . The values y=y + e(0, 1) correspond to the descending 
(transmitted) radiation, while ju - --ascending (reflected) 

radiation. The single sphere S*T -C-l,]jx[o, 2 ij --set of all the directions 

with and ^6Cb,2Jjj $?•[!>» ljx(p,2/j --hemisphere 

for the descending radiation; --hemisphere for the 

ascending radiation. 

In order to simplify writing everywhere, below we will omit the 
vector magnitude sign ("-»■" on the top) of the direction vectors (for 
example 9 T»( 31 */•}•$'• ?'* uishtuiv and the 

coordinates (for example * ,% **,*{* Below we 

will also not constantly mention that — •‘<*»y***. The integrals on the 

single sphere or the hemispheres will be written as: 

In order to reduce the writing in the boundary conditions of the boundary 
value task, we will write 9/o instead of and 9f0 instead 

of With integrated differential equations of transfer in 

the boundary value task, we will omit description of the area of their 
analysis, recalling that everywhere it will be the same: 

«fi* «-0,tf€flT| *4, 4tf.fi *J . 

~~ t 

On the hypothesis of stationary state of the medium and constant 
external radiation flow, the field of the quasimonochromatic polarized 
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radiation is completely described by a four dimensional vector 
Its components in the SP-presentat ion [14] are Stokes parameters 
C * <?«.* which have dimensionality of radiation 

intensity [15, 16, 17]. These parameters are used to reflect all the 
characteristics of radiation polarization, and the first component 

(t , coincides with the energy brightness, or complete intensity 
of radiation at the point in space with coordinate r in the direction s. 

If the medium is macroscopical ly optically isotropic, i.e., there /6 
are no phenomena of linear and circular double beam refraction and di- 
chroism, then the matrix describing the photon scattering law on the 
particle with regard for polarization of the light wave fi(i,&'—a) is 
a function of the scattering angle x [18 - 21]: 


ra,s‘-~s) » /V«, si% fid, e#jcj = fici/wy. - [ty (f,c*x)], 


where 


cax * /*p'+ 'it-y*' co\(v-v) 


We will examine the atmosphere adjoining the horizontally hetero- 
geneous Lambert surface (task for planetary atmosphere with Lambert bot- 
tom). Reflection from the bottom is diffuse, while the reflected light 
is nonpolarized. The matrix of the Lambert reflection is t* -£4YV-i) S(/ -l)j , 
*’,/ -1,2, 3, 4, while the Lambert law for polarized 

radiation can be written as 

( 1 ) 

where , i.e., m £ H,*), 0 , 0, oj t the diffusely 

reflected intensity It* (1* f at each point \ of 

the plane at level is isotropic for all directions , The 

operator of reflection 
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ri - Zl/c^H.soycts'. 


( 2 ) 


Under these conditions complete Stokes vector function 
is solution to the following vector boundary value task 1: (3) - (5). 

$(*) --here and everywhere below is the Dirac function. 

The transfer operator 


fa, ^ ) a * Sin H<UX 


(7) 


The integral of collisions, the function of the source 

« ZscCVJl f><W,y, ft*) +<*,*, ( 8 ) 

&ut&> and --altitude distributions of volumetric coeffic- 

ients of the extinction and scattering respectively. The nucleus of the 
integral operator (8) is the angular matrix onr)f£ 

whose structure and properties are described in fair detail in publica- 
tions [21 - 25]. We note merely that the matrix P in contrast to the 
scattering matrix P is not a function of the scattering angle, but 
depends on three angular variables: >»,/*' and 1 /7 

The albedo of heterogeneous Lambert surface is presented in the form 
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y H Ck) *9>r(%i)*l. 


(9) 


where the parameter -certain constant component which does 

not depend on the coordinates x, y; describes the horizontal hetero- 
geneities , u 

Thus, the task is solved of dissemination of radiation with regard 
for its polarization in a scattering and absorbing layer which is hetero- 
geneous in height with horizontally nonuniform diffusely reflecting bottom 
(without multiplication). 

The boundary value task 1 is linear, and the complete Stokes vector 
can be presented in the form of superposition of two components 

* $?*<*,»)+ (+, 3 ). ^ 1 0 ^ 


The vector p* <?"» tf # * V*J corresponds to the nonscattering 
radiation from a monodi rect iona 1 external source, and satisfies boundary 
value task 2: 


P/ 0 ~T$ a S(3-m6J (11 - 13) 

where the one-dimensional differential transfer operat 0 r 

0* * + {****)> (14) 

and is computed analytically [23]: 

$7?, a)* fj$ A e Tf Q* ) 0,9 €&-]' ( 1 5 J 

The optical thickness at level z , 

a 
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( 16 ) 


The complete optical thickness of the layer z+atffl . 

Component % V corresponding to the Stokes vector 

with regard for repeated scattering in the layer and diffuse reflection 
from the bottom is found as a solution to boundary value task 3: 

?l» m (17 • 19) 


where 




2 . Mathematical Model of Vector Spatial Frequency Characteristics 

Solution to boundary value task 3 is found in the class of general- 
ized functions which assume the existence of an integrated Fourier trans- 
form [ 2 ] in the form of a series of the perturbation theory for the para- 
meter 4 : 

*<• >•> • .***>&.*•«•«* 

(23) 

We substitute (23) and (9) into the boundary value task 3, and pre- 
sent for it the equivalent presentation, boundary value task 4: 


t 


o 

4 


8 



(24 - 26) 


zV»£ - *ilF, 

4r*o *** 

‘0, 

K-O 

£ i” £/„ ~ K , 


from which we will obtain boundary value tasks 
tions of of the solution ^*'in the form of a 


for different approxima- 
series (23). 


With k = 0, we obtain boundary 
boundary values which do not depend 


value task 5 with 
on the coordinate 


coefficients and 
: 




$h u °J ZL 


(27 - 29) 


This is a standard one-dimensional vector transfer 
solution yields a horizontally independent component of 
$£ fe*s)~{T 0 ,Q;V.,V , J » through which illumination 7f£ 


equation whose 
the Stokes vector 
is determined 


£(*)* RT. + 


(30) 


The method, algorithms and programs for solving the boundary value 
task 4 and processing of the calculation results have been well studied, 
tested and presented in detail in publications [22 - 28]. Results of 
numerical solution to this task for a set of optical models of the atmos- 
phere with underlying surface of calm ocean and cloud cover type are pre- 
sented in publications [28 - 30] and reported at the all-union conference 
[31, 32]. The developed technique and the program system AP-5 [25] make 
it possible to calculate the Stokes vector with sufficiently high accuracy. 
With physically correct optical parameters of the atmospheric model, rapid 
convergence of the iterations and stability of the finite differential 
method are observed [28]. 

We see further development of the technique for solving 1 the task of 
dissemination of radiation in planar media with regard for polarization 
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in the need for developing a method for calculating the diffused reflect- 
ion from uniform and nonuniform underlying surface based on solution to 
a one-dimensional vector transfer equation for an insulated layer. This 
is possible with generalizing the SFC method [5, 6, 33, 34] for the vector 
task, i.e., with creation of a mathematical model of vector SFC. 

We will introduce a designation for the Fourier transform operation /9 
for the coordinate : 

&CW) m JLf rt4) * rAt <^ > *4 (31) 

The f function can also have other arguments which differ from \ . 

The spatial frequency P-ih'tyj only adopts real values. The inverse 
Fourier transform is defined as: 

In order to reduce the writing, we will indicate everywhere one sign of 
the integral, recalling however, that integration is done by pairs of 
variables (x, y) or(P x » p y). We will introduce the spatial Fourier spectra: 
for albedo variations: 


1( f )* TC<b3<P> > (32) 

for the Stokes vector component (we designate the Stokes vector components 
by printed letters, and the Fourier patterns by the corresponding hand 
written letters): 

jy K djP*, 3 ) 1 &~CZJ(fk); G (•,/*,•’) * STS *HM; 

' Z'CVJ(P')-, V(t,Wh Sfv.3 (V, ( 33 - 36 ) 

For the Stokes vector 

i 

(37 - 38) \ 
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k*I --index of the expansion term (23). 


The Fourier transform (31) applied to expression (6), results in a 
linear operator, which is one-dimensional in space: 

j((f>) s)= ( 39 ) 

In this case we use the ratios [2]: 

' *,‘>e **?>_& [ft '*'$«> 


The integrated operators of the function of source $ and reflection j? do 
not depend on the coordinates ^ , therefore 

v](p) } chi (p). 


With k = 1, from the boundary value task 4, we obtain boundary value 
task 6 for linear approximation of a number of perturbations (23): 

'$$=3%} file* 0 ) + (42 - 44) 


After Fourier transform (31), in boundary value task 4 we arrive at 
boundary value task 7: 


~K/.=0> 


*«**>*' 


(45 - 47) 


in which wave number Pj is a parameter. We will attribute to wave numbers / 1 0 
P, the index corresponding to the order of approximation in expansion (23). 

We note at once that all the boundary value tasks with differential oper- 
atPr f((p) are vector transfer equations with complex coefficient of 
extinction equal to t(i and consequently, the unknown functions 

will be complex. The numerical method of solving the transfer equation 
with complex function was developed [36]. Its generalization for the 
vector equation does not represent a problem, since it naturally is com- 
bined with the method developed for solving the standard vector transfer 

o 

4 
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equation [22 - 25, 28]. 


It is easy to show (for example, as done in publication [5]), that 
because of the linear nature of boundary value task 7, the Fourier pat- 
terns of vector and albedo are linked by a linear relationship 




(48) 


1 [Xid*,*) 1 

/= / I'EfW, 
Ui a, fa *> I I 5> J 

*i/l (*, ft, 3 V L Tl 


(49) 


where the vector linear SFC 


V* (**.*)• [h,X», Tt ) 


satisfies boundary value task 8 


Kwfi ftl^tfi****) 1 * 


(50) 


The vector Vf does not depend on the conditions of illumination of the 
layer and variations in albedo 9r ' v . The inverse Fourier transform 
yields : 




(54) 


For the linear approximations of the series of perturbations (23) on the 
order«*2 , we obtain boundary value task 9 from boundary value task 4: 

, < 55 • 57 > , 

0 

i 
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in which two successive approximations are recurrently linked through 
the first component 1^ of the vector ft which is included in the boundary 
condition. For nonlinear approximations, the following occurs 
Theorem 


If the vector function 

^ (^) A'y-i /V/ 1 ^ J ' 


(58) 


with fixed values of spatial frequencies Pi,...,P|< satisfies the set of / II 
parametric boundary value tasks 10 which are one dimensional in space: 




then the Fourier-patterns , the vector of the functions 

9 ,s) > being solutions to the corresponding boundary value task 
9, are unequivocally defined through the vector f unction yf* with the help 
of the integral ratios: 




Proof . 


We will make the Fourier transform (31) in boundary value task 9, and 
will arrive at boundary value task 11: 

k(M 

$1=0) = • •' (63 - 65) 


We will make an induction assumption, by introducing vector function 
so that the following ratio occurs 
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( 66 ) 


then 


Ww-I «,*,*)e <<s 34 

\TlftVe <* *>4* •&£?„ *;/fc*,j. 

f /•». ' 


1 «>«*> Xr- t a,^k 1 

• ctp*-4 -. <t Pi . 


*.»> 

ifaat-. 


(67) 


We substitute (67) into (65): 


K(M$! (*>*'*) = &) 3 /,•<); ‘ %/„•(} fo*<fid+ 

« (68 - 70) 

■f 7T~j7 7iw? 

(J) "° 

Assume that the vector function - tyj* ( 2 ft ,3) with fixed para- 
meters Pj P k satisfies the boundary value tasks 10. We multiply the /J_2 

left and right sides of (59 - 61) by the product i (P*-P*-*) 

and we integrate for parameters P^ Pj in limits from -**to + •* As 

a result, taking into consideration that operators 'ftW' , $ and $ can be 
removed as signs of integrals for parameters, we obtain: 


1 K (p*) a 3 J* 

- 4>J* 3 ' /y_ fCf " 

* ^ C ^9 fa-*) fa &) •• ^ ^ 

« ^ ^ 
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-JUa 
. Jp„~H - 

. Xk (H,e*,~>fo M*' ~ l/p J + 


By comparing (71 - 73) and (68 - 70) we establish that assumption (66) is 
true, and occurs for all n*2, while for the vector functions^, presen- 
tation (66) occurs through the corresponding vector function T?'k- Thus, 
if the vector function satisfies the boundary value tasks 10, then the 
corresponding vector function^* is determined from formulas (62). 


We now assume that for all k' = 2,3, ...k, ratio (62) occurs. Then, 
by substituting in (63 - 65) instead of expression (62), and also / 1 3 

(67), we arrive at an equivalent boundary value task (71 - 73). With 
fixed values of parameters from this task we obtain the^ara- 

metric set of boundary value tasks 10 for the vector function |< , i.e., 
if presentation (62) occurs, then the vector functions iff* sati sfy boundary 
value tasks 10. 


Solution to boundary value tasks 10 is determined only by regular 
parameters of the layer and backing model: ♦ 4atto» + 

Consequently, the vector functions are universal characteristics of 
the nonlinear transfer system of radiation with regular parameters, do not 
depend on the conditions of illumination and are invariant in relation to 
the horizontal heterogeneities of the albedo The vector functions 

fa**) are called vector SFC of the k-th order of approximation 
for parameter £ , corresponding to the term *'«» the expansion (23). In 
summarizing what has been said, we obtain a calculated formula 


[ fa'~' fa *) e m ■ ' (74) 
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CZJ) 


ft, 


E . jJq a, p*,. *) <tp * • •• ^ e 


it /.- AW r , 


•jjh******! ■ •&£*«*■■ *«•>*-*• 


Thus, if there is a known solution to the one-dimensional boundary 
value task 10, then any term in the perturbation series (23) can be 
computed through the vector functions fJ* using formulas (74). 


We will return to boundary value task 10. In this task, the depend- 
ence on parameters Pi»***>Pk-l is introduced through the term in the 
boundary condition (H,/w fk »•)’• Since all the operators £ ,3 

and ?in boundary value task 10 do not act on parameters P j , — .Pk-1» then 
because of the linearity of the transfer equations, factorization by para- 
meters occurs: 

$?« (*/ P*r"t Pi* 9 > m ** ( 9, P* (75 ) 


and consequently, (see [58]) 

X K pK, f/, i ) m % L (*/P‘» s ) 8 ) . 


(76)' /U 


By substituting (75) and (76) in boundary value task 10, we arrive 
at boundary value task 12 with one parameter P k : 


%l s0 > 


(77 - 79) 


But this boundary value task is adequate to boundary value task 8: in 

both tasks with the same integral 3 and differential *>£(p)i operators, 
the sources, coefficients and boundary conditions coincide. It follows 
from here that for all **2i- ' 

3 

i 
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(80 - 81) 




with k = 2 


(82) 




(83) 


Similarly, after further continuing the recurrent link with the help 
of ratios (75 - 76), and (80 - 81) and 


RX«- (H, p*v-vM> I fatWfa*), 


(84) 


we obtain a presentation of the vector SFC of any order of approximation 
k>2 to the vector SFC of the linear approximation: 

V K ^5 dfaQfftfo&hA (85) 


This is an extremely important ratio, since it reduces solution to the 
recurrently linked boundary value tasks 10 with increasing number of 
parameters Pj,...,P k to a solution of one boundary value task 8 with one 
real parameter P. By substituting (85) into (74), we find an expression 
for the expansion term (23) through the linear SFC: 


^ 


( 86 ) 



As a result, the Stokes vector can be computed by summing the series / 1 5 
(23), assuming parameter £ -I, in the form of superposition of three 
terms : 

4>(*,S)= %(*,*)+%"’*)+ <?*<*'*>> (87 ) 


where fidt,*,*) --linear approximation of (54); the component of nonlinear 
approximations, starting from second and further: 




- -solution to the one-dimensional vector boundary value task 5; 
the quantity E (30) is determined through the first component Iq(z,S) of 
vector i f (P)--Fourier pattern (32) of the variable component of 

the bottom albedo l the vector 4*£^>A*jJ--solution to the one-dimen- 

sional vector boundary value task 8 with .parameter P k ; the quantity 

Cx(Pe) s KXi(H,Pt,s) is calculated by integration for the angular vari- 

able S (2) of the first component Xj of vector^ only at the bottom level 
(z = H) on the hemisphere of directions •%£&•' 


3. Vector Functions of Scattering, or Vector Functions of Influence 


We will examine the series of perturbations (23). We will introduce 
the vector function of scattering [3, 4, 6] of the first order .fcoejXu,*) 
so that 




(88) 


i.e., the vector scattering function is the inverse Fourier pattern of 
the vector linear SFC ip £ 



&(*/***> m (jnpJEl *4 (*>?>*)& 


(89) 
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With the help of ratios [2] 



( 90 ) 


(91) 


By inverse Fourier transform in boundary value task 8, we obtain boundary 
value task 13 for the vector scattering function of the linear approxi- 
mation 


$!(*)**, s) 0m, On) i 


(92) 


fou)) 5 • 


(93 - 95) 


By comparing this task with boundary value task 6, we are convinced that 
is the Green vector function of boundary value task 6 for 
linear approximation of the Stokes vector . Or in terms of 

the monograph of V. S. Vladimirov [2], a) is a fundamental solu- 

tion, or a vector function of influence of the vector boundary value task 
6 with perturbed boundary condition, more accurately, with isotropic^ - 
source on the boundary Z = H at the point with coordinate 


/ 16 


In order to conduct the numerical calculations, we give preference 
to the following algorithm: initially boundary value task 6 which is one- 
dimensional for space is numerically solved by the vector function 

with a set of parameter values P; through it using formula (89) by the 
method of rapid Fourier transform, the vector function of scattering 

is calculated. Numerical solution to boundary valye task 13 
directly for vector function is difficult to realise on a 

computer, since this task contains five variables (w^r) 



However, in order to solve task of the boundary value type 13 containing 
local spot or spatially limited sources, method of SFC is also applic- 
able which still has an advantage: in these tasks, the series of the 

perturbation theory of type (23) only contains linear terms. An indi- 
vidual publication will cover this problem. 


The Stokes vector ^ (54) with the help of (88) can be presented in 
the form of an integral packet: 

; f(#)e . 


• uD <hn e~ £( fc H'/v • C '<*+*, *)* 

<[c&l* l( *' ***'*+]+*» * * 


(96) 


We will use replacement of variables ‘ 

formula [2] 

We note in particular, that the following ratio follows from (88): 

^£*(2,3)2 Vi(*tP a0 t 8 ) (98 - 99) 

from which follows the existence and the limited nature of the integral 
in the right side of (98) because of the limited nature of the vector 
function • The vector function is real and is solu- 

tion to the one-dimensional vector boundary value task 14: 

4 $?*« 35 ?/; m.'°> (100 - 102 ) 

1 


, and also / Y7_ 

(97) 
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Evidently 


2*i <H*A*) =£#e Kl (H, , 


(103) 


R9 Al (i^H,s) (H, r,s)e'‘ ( Mdf>. 


(104) 


We introduce the vector functions of the k-th order of approximation 
of series (23): 

Gk ( i, Tik,..., Zii, o) (*> far >fa tfe'* (105) 


or 


lf K (*/ far ■■> fat a J^Gtc 


and we substitute (105) into (74): 


(106) 


i,z) a j£j^}Jl(fa)Hfa-Pihf(fa-fa-d e # > 
QJi,tl2K,."> W,»; ' Ifl****" 

mB ff % (U*>e (h,%i hi fr (tu *> e (fk $>(**)• * 

*»«,*)**&& 

-l (fal*)] dpi dtitK ••• Ail * Ejh/fyCkl)— Cttu) * 

(J%g) 

.&(*, «.»/,«;• [JarTyrR ‘fh ***** *-*0*]’ 

* [cirfr fa'* (fa/lu'tti +b*+j)] dpt]* 


-x iK +iii~)idp^y[fajjxfC&&-+++»w* • 

« Aa — Atar « 


A' 


( 107 ) 



~E J"J { r k*)(^XlK"'^XU ^ 

~w (,t> 

* l(T x /-t4z + tjjf*/)$(tu-'k*+ %K*Jt)"' iCkt-i-TkK* %»k-i)* 

* 8 * ‘k + £«) cflx2K "* frlK+I * £ y ^/ '%(***)•• (It*)* 

*(&(*, ^ -\* , Tt* - W ^ ' *<> - «**« * 

= £ jJtf (Xtz- Xu) fjjfrs- Xu)" 9r <%-%«)• 

— 

* 5c f*> *«/ •» Xu, 5) 

It follows from (106) 

^ (g, 0,-, q *)«//§,(*, 


/IS 

(108) 


Using presentation (85) for the vector SFC and formulas (89), (103), 
and (104) for expression (105) we find the most convenient form, factor- 
ized, split for parameters P^: *.< 

$ (i ’ T “ > K " ’ s) "(h)“j3 ¥ * c *' **< 

texf >[~i £ (/>*, <*? i * /~zfev£ H <*/A/ ^ ^ 

irhi Lh L (»,?‘,*) z tiF ‘ f %]= 


(109) 


« fl [w 
1=1 




- (\*,*>*) J] f $&** Ckt, //,&). 

As a result, the vector scattering function {?* of any order of approxi- 
mation k» 2 ' of the perturbation series (23) can be defined through the 
vector • scatter ing function of linear approximation ff2 It is further 
easy to show that any term in the series (23) can be found using the 
integral packet: 


£ (Xt t *, 4>* 9Ak*~ 

k) $lCUk, t, 0^1 $0n(1te,Nr3)dhA ...</♦ H ■* 

,E J %(ki)" fy/ 1 **) & <*-*•>*>*)' 

^ (ti *‘ “ k *> H > * v<k< - v * 1 

' . i 

0 

4 
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The Stokes vector can be computed by summing the series (23) similarly / 19 
to (87), in which the linear approximation*^* is computed according to 
formula (96), while 

/ fl / dll) 


The vector functions <9* satisfy boundary value tasks 15, which are ob- 
tained after an inverse Fourier transform in boundary value task 10 using 
the relationships (90 - 91): 


(rfe* («. v, »). sr>»)«e^> 


(Ilia - 111c) 


where the differential operator 


in.) (*>%*)**•"* < llld > 

acts on the variable , while the variables in this task 

are parameters. We establish from here that factorization occurs coin- 
ciding with (109): 


R0 tlK ft C*,kK, *)?[}$« (Xtc, H, *) . 




4. Calculation of the Constant Component of the Albedo in Vector SFC and 
Scattering Functions 

We return to boundary value tasks for the vector SFC . The multi- 
plier in the boundary condition (53) only depends on the parameter P. We 
will designate it 

t 

tl# = ~ ‘j j 12) 
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Because of the linearity of boundary value task 8, the vector function 


•> = £ w *> ^ 3 r 


introduced with the help of the ratio 

C*,h *) * t CfO W (*,**), 


will satisfy boundary value task 16: 

fif) $ 1 ,- 0 ) *iw. 


(113) 


(114) 


(115 - 117) 


After a chain of several transforms: 

l(f>)* ft if) 7 (Hfh*) (118 - 121) 


we arrive at expression 




SSiihJL, 


( 122 ) 


as a result, the vector SFC of the linear system of transfer of radia- 
tion "atmosphere-uniform Lambert underlying surface with albedo q" is / 20 
analytically expressed through the linear SFC of the isolated atmos- 
phere. We note in this case that 


c<fi* 


(123) 


is determined through the first vector component *7, and only depends on 
the spatial frequency P. Now the vector SFC (85) of any order can 
be computed according to the formula 




ff* e(M 


( 124 ) 



Factorization for parameters is a convenient circumstance for realizing 
numerical calculations reducing the volume of computations. The compo- 
nents of the Stokes vector (87) can be computed in this case by merely 
knowing the vector SFC of the linear approximation W for the isolated 
atmosphere: . 

■***&&%$ §£rk, 

r, • . . .. ( 125 - 126) 


ft)- P*-i)exr&* (/*,*>]< 

eUf 


It is apparent from these formulas that dependence of a Stokes vector 
on the constant component of albedo is not linear. However, it is easy 
to obtain from these same formulas an evaluation of the deviation from 
the linear relationship, by distributing the fractions l/l-qc(P) into a 
Taylor series according to the degrees of argument qc(P). Using (122) 
for the vector scattering function one can also obtain a formula of 
linkage to the constant albedo component: 

(127) 


©(«*>•-> ‘e&J. 


After this, the expression for can be written as 


(128) 


However, in this case, we do not succeed in presenting the Stokes vector 
in the form of an integral packet similar to (107) which would include 
only the scattering function determined for the isolated layer. 


We introduce the vector scattering function &* s [9**i9y**9**>&r*J 
through the SFC vi£* , taken with zero value of q: 


Starting from (122) and (124), one can write 


25 



= fot ft, g ; ^ a) ^ 

<P* ' ' V- C( ft )) y . , St. $) = c(f)t)' 


For the vector scattering function y,' - j 0a ; $nr / $ **,&*} we have 

& ( i, st « .Az + * ft *», a) e*t £■'(/ %* *^St, 

e7W,/ 

We use expression [74] 


/ 2 1 


^ ^ 2 , s) = j &*-/>-') e 




r K ' +i~i-f (ay* 1 

ft* s) 


1 ra-w* dr ‘~ dri ° ««“ 


n WMfohyfcH 

Sj&HZ ‘ C^n^e. -A/ 







'A 


r * * 


JSSU-z-c 


This is not reduced to an expression of the packet, since the integrals 
for the parameters contain functions which differ from the constant, and 
therefore they are not taken analytically and do not result in ^-functions , 
as occurred, for example, in the case of derivation of formula (107). 
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Here the conclusion is suggested that a sequence of calculations be orga- 
nized in the case of nonzero component of albedo q (if q = 0, then the 
just derived formula is reduced to formula (107)). Initially we calculate 
the vector linear SFC as a solution to the boundary value task 16. 

Through it we compute the vector SFCijf^ according to formula (122) and 
the linear vector scattering function according to formula (127), while 
the corresponding nonlinear approximations are defined from formulas (124) 
or (85) for ^> , and according to formulas (128) or (109) for & * 

5. Calculation of Diffuse Reflection from a Uniform Lambert Surface in 
Solving the Vector Transfer Equation 

From the concepts obtained above for solving the vector boundary 
value task 1 with horizontally nonuniform boundary condition, described / 22 
by the Lambert reflection law (1), one particular, but extremely impor- 
tant result follows. We will examine boundary value task 5. It differs 
from boundary value task 3 because in the latter, the albedo is horizont- 
ally nonuniform. We did not place any restrictions in presenting albedo 
W in the form of (9), except one: there should be an integral Fourier 

transform (31) for the component %(Xt) . We place in (9) q = 0, while 

-<J-const . in this case the Fourier pattern 

= TC<tnl(P) =_X. <tn expCHfi - (&)*■ %(?), (129) 

The linear vector SFC 

^ ft = w (130) 

The linear approximation of the Stokes vector (125) looks like 

The formula of nonlinear approximations of Stokes vector (86) is simpli- 
fied: 

i 

.i 
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£(?,s)= E J'-f S(Pt)8(fr -Pi)" 

tewfcfaiflfy-JwfyefaNtor+fr *♦; fif'S (*,*i 


( 132 ) 


The vector function 


iv,(*» s ) a wC»;p=A^*f^i^/ t/ »*^ J (133) 

satisfies the standard one-dimensional vector boundary value task 17 with 
isotropic source of nonpolarized radiation on the lower boundary: 

ifr-lV.) . ( 134 - 136) 

According to the principle of the maximum (35) and the properties of the 
operator'.^, the value 


C, a C(f>*6) * C W ,*) H- 

We will sum series (23) in the form: 


(137) 


« ?/(*, a)* a, s\,' ( 1 38 ) 

^ vij --solution to the boundary value task 5 with albedo q = 


0, 


? r £ «»*> - f »<• 


(139) 


The series for ^ was summed as a geometric progression with denom- 
inator , since usually f/< I,. The value of the magnitude of 

illumination because of diffuse reflection from the uniform Lambert 

I 

surface can be obtained by another method, based on presentation of terms 
in the series (23) through the vector scattering function. With %Ck) 

= q = const, yje find from (96) with the involvement of (98): 
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fy*,a)= fta, *, *£>#,/>«* a >^ V/(t,8h U 40 ) 

from (107) and ratio 


(*> 4 -m ° > &,£) &&> *v~» H$V*« ~ Au (141) 

following from (106), we find that 

- etS £l& a, xu'rtiu, *)'**- *tu ■ * 


But according to (85): 


/23 


Vk(2/Qz ~ (*»4)^iJ f 


(143) 


where 

value 

meter 


the real vector function z/V satisfies the standard vector boundary 
task 18 obtained from boundary value task 8 with value of the para- 
P = 0: 

T®*/.*®'’ (144-14 


while the quantity 




(147) 


We will sum the series (23): 


*-fr<3r # 


(148) 


The boundary value task 18 is a particular case of boundary value task 8 
with P = 0. Consequently, for vector ffi* , the same properties are cor- 
rect as for $j£'. In particular, similarly to (114) we have: 
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(149 - 151 




Thus , 


** i- *rX# 


(152) 


Assuming q = 0 (compare (152) and (139)), we obtain 

(153) 

It is appropriate to note here that the obtained formulas for calculating 
the horizontally uniform albedo component result in two important conclu- 
sions. First of all, for. calculations of the Stokes vector of the 
atmosphere for the uniform Lambert surface, it is sufficient to solve the 
vector boundary value task 5 for the vector <f% G,8) with zero boundary 
conditions (in (29), we assume q = 0), i.e., for an isolated atmosphere, 
and then we take into consideration the component of the Stokes vector 
^ (139) governed by illumination because of diffuse reflection from the 
bottom. In this case, the quantity JE is computed through the first compo 
nent of Stokes vector and parameters of external parallel stream on the 
upper boundary (30), and determine the illumination of the atmosphere 
bottom. The vector function does not depend either on the illum- 

ination conditions, or on the albedo of the bottom, as a universal charac- 
teristic of the atmosphere itself, essentially describing the propagation 
in the atmosphere with regard for polarization of diffuse natural non- 
polarized light emitted by the uniform backing. Under conditions of obser 
vation from above, the functioning practical ly coincides with the charac- 
teristics of atmospheric- transmission with regard for the effects of 
polarization and multiple scattering. The coefficient Cq is an analog of 
the spherical albedo of nonuniform atmosphere illuminated from below by 
parallel beams (9). It is obtained as a result of integration of the 
descending radiation of all directions of dissemination , 

and in all directions of incidence of the external stream from below 

4*€5T on the boundary Z = H. Secondly, from the obtained results. 
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our hypothesis regarding the fairly random solution of the constant compo- /24_ 

m 

nent q in presentation of the albedo in form (9) is confirmed: in addi- 

tion to the quantity q, a certain constant component can also be con- 
tained in what we call the albedo variation, i.e., in the quantity q, r (Oj . 

Its contribution is taken into consideration analytically using the 
component of Stokes vector ^ according to formula (152). 

6. Evaluations of Nonlinear Approximations 

* 

Using the presentations (96) and (107) of the terms in series (23) 
through integral packets of albedo variations with vector scattering 
functions, it is easy to obtain evaluations of the contribution of non- 
linear terms with numbers ,K 4 *?.,--any assigned number. If we assume 

0* 9^ • max 19,(101 ^ .1 (154) 


then from (107) with regard for (105), (85) and (143) we have: 

... ct%a ✓/v... ys K o , a) • 

. - £<*% Q (*;s) (2x/(n,s*) K -f 


Evaluation from above of the contribution of nonline*’" approximation (23) 


s 

’ t £ & (for*- 

•Jfc/ ^ **** £ W.C*.*)***^ 


(156) 


With q = 0, the formula is simplified 

4 *, t 6 <& *"' 4 * 


(157) 


i It is important that evaluations of nonlinear approximations are 

expressed through the characteristics of an isolated atmosphere, which can 

% 
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be calculated using the simplest boundary value tasks which have already 
been discussed above. With k^ = 1, we obtain from (156) the upper eval- 
uation for complete solution to the vector boundary value task 3: 






1- ( ++ $•***) c» * 


(158) 


where --solution to boundary value task 5 with zero boundary 

conditions. 


7. Fundamental Solution to the Vector of Transfer Equation with Local 
Perturbation in the Boundary Condition 


If we take *5^ , then f(P) = 1. We find the expression 
for the corresponding Stokes vector ft, based on the presentation (74) 
and determination of (106): 


♦ ®K <p,,Xte)]4$ K . 

oo 

- (*/ &*>•••/ SfXt* - ... « 

= E 9 k (*/ % ?>-> 


(159) 


/25 


Here the following designations have been introduced: 


C&)*L 



-a*r- 


c x fye)* fan 


(160 - 161) 


Based on presentation of the vector ft in the form (110), after simple 
transforms, we arrive at the same result: 

_ ft CM,* 

-iOtodbC*'**)*** ( 162 ) 


\ 
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It follows from this that the vector scattering function 
and \ 0 r'j^s) are fundamental solutions to boundary value tasks 

6 and 9 respectively with S -source on the lower boundary. As a result 
of summing series (23), we obtain an expression for the fundamental solu- 
tion to boundary value task 3: 

•&<*> *>♦ £ ^ / 

&,$) + ) eq» m cc-p)~ C{m0, (163) 


8. Amplitude and Phase Characteristics of the Vector SFC 

It follows from what has been said above that in order to solve the 
vector boundary task 3 with regard for any order of nonlinear approxi- 
mations of the perturbation series (23), it is sufficient to solve the 
task for the vector SFC of linear approximation ii for an isolated 
layer, i.e., initially solve the parametric boundary value task 16, and 
then calculate according to formula (122). We introduce the amplitude 
and phase characteristics through the diffuse component Jf of the vector 
SFC 







> l <r& 





m *1,1,1, '*■ 


(164) 


The vector AFC is S)s[Am ) f --vector FFC: 

Am = i| (te Xm)S ' } ‘P/n * anctgOfmXm/RtXm] / Xm*0. (165) 

It is easy to establish that 

with observation in the nadir of no phase distortions. m ' 

\ * * / 

vector of shift in tjie horizontal plane with inclined sighting routes: 

1 
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--Bouguer attenuation. The vector 




which / 26 


is determined for those components with number m for which A m ( 0 ) ^ 0 
coincides with the frequency-contrast characteristics FCC, which only 
takes into consideration the amplitude distortions, A° m (P = 0) = 1. The 
introduction of nonlinear AFC and FFC is surplus, since all the nonlinear 
SFC^are expressed through linear SFC vfc. The introduction of AFC and 
FFC for the- vector ‘^T i s less convenient, since the calculated formulas 
become less graphic, since c(p) is not always a real function. 


Sometimes it is convenient to introduce the scattering function 
thrugh the diffuse part of the linear SFC: 


f f*, *,*,*). 

The complete linear scattering function then is 

' A ‘ fCk ' X>t " " 

9. Direct Transfer Operator of the Atmosphere with Regard for Radiation 
Polarization 

In summarizing, based on the results obtained above, one can write 
an expression in which the Stokes vector of radiation reflected by the 
atmosphere-underlying surface system will be analytically linked to the 
albedo 

Cp = ff (<bO~Qo («,«) ♦ A<\ '/#) ♦ vjr 

(170) 

--solution td boundary value task 5 with zero condi- 
tions. Illumination because of constant albedo component q (139) is 


(167) 


(168) 


(169) 
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computed through the quantity for isolated atmosphere: 


¥♦<*/*>* (171 - 172) 

--illumination of the bottom with incidence of parallel rays on the upper 
boundary in direction So; --spherical albedo; Wq(Z,S) = 

--vector transmission of isolated atmosphere illuminated 
by parallel rays from below. Illumination ft (30) with q ^ 0: 

E= (173) 

The linear approximation for calculating the horizontal heterogeneities 
of the albedo (54): ^ . 

(174) 

The contribution of nonlinear approximation is computed / 27 

according to formula (111) through the linear vector scattering function 
or according to formula (87a) through the vector linear SFC $£ Upper 
estimates of the contribution of nonlinear approximations are obtained 
from formula (156). The direct transfer optical operator of the atmos- 
phere (170 is a generalization of the similar operator for the 
scalar task [13]. 

The authors are sincerely grateful to M. V. Maslennikov and A. A. 
Ioltukhovskiy for useful discussions of the results of the work. 
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